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Abstract 

We study the gluon propagator in Landau gauge in the deconfmed phase of 
SU(2) gauge theory. From the long-distance behaviour of correlation functions of 
temporal and spatial components of the gauge fields we extract electric (m e ) and 
magnetic (m m ) screening masses. For temperatures larger than twice T c we find 
no additional temperature dependence in m e (T)/T, while m m (T)/T drops with 
increasing temperature. The decrease is consistent with the expected behaviour, 
m m {T) ~ g 2 (T)T. We find m e (T) = 2.484(52)T and m m (T) = 0.466(15)# 2 (T)T. 



A basic non-perturbative feature of the high temperature plasma phase of QCD 
is the occurrence of electric and magnetic screening masses for the gluon. They 
play an important role in controlling the infrared behaviour of QCD Q. While 
the electric screening mass can be calculated in leading order perturbation theory 
and is found to be O(gT), little is known about the magnitude of the magnetic 
screening mass, which is expected to be 0(g 2 T) |]J. However, also other functional 
dependences have been obtained in approximate non-perturbative approaches [§, |3| . 
Recently some attempts have been made to determine the magnetic screening mass 
through the analysis of gap equations [|], |[ ||. 

It is the purpose of this paper to discuss a detailed study of the temperature 
dependence of the magnetic screening mass, m m (T), obtained from calculations of 
the finite temperature gluon propagator on the lattice. Early lattice studies have 
attempted to extract this quantity from gauge invariant observables, using lattices 
with twisted boundary conditions |7], |J. Likewise the electric screening mass usually 
is extracted from gauge invariant Polyakov loop correlation functions. Here we will 
determine the effective gluon masses directly from the long-distance behaviour of 
the gluon propagator. Although this requires the fixing of a gauge and leads to 
a gauge dependent definition of a screening mass, it has the advantage that it is 
closest to the perturbative definition of these masses and allows a direct comparison 
with perturbative calculations. A more detailed analysis of the momentum space 
representation of the gluon propagator may, however, also lead to the determination 
of a gauge independent pole mass 0. 

We have performed calculations on large lattices of size N% x N T . In most cases 
we use N T = 8, which insures that the calculations in the high temperature phase are 
performed at gauge couplings well inside the scaling region of the SU (2) lattice gauge 
theory, i.e. (3 > 2.6. For the spatial lattice we generally use N a = 32, which allows 
us to analyze correlation functions up to distances zT = 2. We have performed 
calculations in a large temperature interval from T ~ 1.3T C up to T ~ 16T C in 
order to become sensitive to possible logarithmic corrections to the leading linear 
dependence of the high temperature screening masses on T. Typically we have 
generated 800 independent gauge field configurations on which we fix the Landau 
gauge, Ic^A^x)! 2 = 0, in order to analyze gluon correlation functions. For the 



gauge fixing we use a combination of overrelaxation and FFT algorithms [QQ, 11 . 
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We have discussed the performance of this algorithm in Ref. fT2||. 



We define the gauge fields, ^(x), in the usual way from the SU(2) link variables, 
[/^(x), the dynamical variables in the lattice formulation, 

A,(x) = g (^(x) - C/t(x)) . (1) 

In a perturbative context the electric and magnetic screening masses are defined 
through the zero momentum limit of the gluon polarization tensor, IL^Qoo, P), at 
vanishing Matsubara frequency po- 111 coordinate space this limit can be realized 
through the long distance behaviour of correlation functions of gauge field operators 
in one of the spatial directions of the lattice, e.g. in the ^-direction, 

G,(z) = (TrA fl (z)A,(0)) , (2) 

where the fields A^(z) = Y^x , Xl ,x 2 ^(^o, are obtained by averaging over a hyper- 
plane transverse to the z-direction. In this way we project onto po = p\ = pi = 
and the limit p% — > corresponds to z — » oo. The long-distance behaviour of 
G e (z) = Gq(z) ~ exp{— rh e z} is related to the electric screening mass, while 
G m (z) = 0.5(Gi(z) + G2(z)) ~ exp{— m m z} yields the magnetic massQ. Due to 



the Landau gauge condition, it is easy to see that Gs(z) = const. JT3 



A problem which may arise with the fixing of the Landau gauge is the occurrence 
of the so-called Gribov ambiguity flilfl . The Landau gauge condition |9 M yl M (x)| 2 = 



is realized on each lattice configuration by maximizing the trace of the link fields, 



£ = [Up(x) + Ul(x)] , (3) 



in the space of gauge equivalent fields |jl0| , |TT], I2|. This maximization procedure, 
however, is not unique. If one performs on an ungauged configuration first a random 
gauge transformation and then maximizes E, it is possible that one ends in a different 



a We denote by fh masses in lattice units. As all calculations, except one, have been performed 
on lattices with N T — 8 sides in the temporal direction, the masses in units of the temperature are 
given by m/T = N T m = 8m. 
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Gribov sector. To investigate the dependence of m e (T) and m m (T) on Gribov copies, 
we took 100 configurations, created on each of them 25 random gauge copies, and 
then performed the gauge fixing. Indeed, we find different Gribov copies, which 
are distinguished by the value of E. In order to test the sensitivity of m e (T) and 
m m (T) on this we selected the two samples of 100 configurations with smallest and 
largest value of S, respectively. On each sample we calculated the screening masses 
from the exponential decay of the correlation functions G e and G m . Within our 
statistical accuracy we do not see any dependence of the screening masses on the 
different Gribov sectors. 
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Figure 1: The electric, G e (z), and magnetic, G m (z), correlation functions for various 
values of the gauge coupling as a function of z calculated on lattices of size 32 3 x 8. 
The temperatures corresponding to the different f3- values can be found in Table 1. 
They cover the interval T/T c e [1.5, 16]. The data points in the upper figure have 
been displaced horizontally for better viewing. 



In Fig. [I] we show the electric and magnetic gluon propagators at various tem- 
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peratures. Already from this figure it is obvious that the two correlation functions 
lead to quite different temperature dependences of m m /T = 8m m and m e /T = 8m e . 
While m e /T does seem to be temperature independent, the magnetic mass clearly 
rises slower than linear in T, i.e. the correlation function G m (z) becomes flatter 
with increasing temperature, which suggests a decrease of m m /T. 

We have analyzed the long-distance behaviour of the correlation functions by 
studying the behaviour of local masses, 



Gi{z) cosh (rhi (z - ^ 



G t (z + 1) cosh (rhi(z + l-^ 



e, m 



(4) 



as well as by fitting the correlation functions with the ansatz, Gi(z) = c cosh(mj(,2 — 
N a /2)). From an analysis of the local masses we find that they slowly approach a 
plateau for zT > 1. This behaviour is shown in Fig. for (3 = 3.12. 
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Figure 2: The local electric and magnetic screening masses for (3 = 3.12 in the range 
zT e [0.125, 1.875], i.e. 1 < z < 15. 



Also the single-cosh fits yield acceptable x 2 " vames om y) if the short distance 
part of the correlation functions (zT < 1) is left out from the fit. This is obvious 
from Fig. [l|. We generally observe that a fit to the propagator for zT > 1 yields 
masses which are consistent within statistical errors with the local masses extracted 
at zT = 1 (z = 8). In Table |T] we give the results from our fits at all values of 
the gauge coupling analyzed. These results are also shown in Fig. || In order to 
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p 


T/T c 


m e (T) 
T 


m m (T) 
T 


2.60 


1.32 


2.28(15) 


2.42(10) 


2.74« 


1.33 


2.83(27) 


1.91(12) 


2.65 


1.53 


2.25(17) 


2.25(9) 


2.74W 


2.00 


2.69(10) 


1.92(3) 


2.88 


3.03 


2.36(14) 


1.63(5) 


2.97 


3.93 


2.26(15) 


1.46(5) 


3.12 


6.01 


2.90(22) 


1.24(4) 


3.20 


7.53 


2.40(15) 


1.17(3) 


3.34 


11.10 


2.26(14) 


1.10(4) 


3.47 


15.88 


2.57(15) 


1.06(4) 



Table 1: Electric and magnetic screening masses in units of the temperature, 
m e (T)/T and m m (T)/T for temperatures 1.32 < T/T c < 15.88. Calculations have 
been performed on lattices of size 32 3 x 8 except for the cases (1) and (2), where 
calculations have been performed on a 32 3 x 12 (1) and a 32 2 x 64 x 8 (2) lattice. 
The masses have been obtained from fits to the correlation functions for distances 
zT > 1. 



relate the gauge coupling used in the calculation to a temperature we make use of 
a parameterization of asymptotic scaling violations of the SU(2) /3-function [13 . 



We clearly see that the electric screening mass is proportional to the temper- 
ature. Up to temperatures of 16T C we do not see any indication for the expected 
perturbative behaviour m e ~ g(T)T. Averaging over the results for T > 2T C we 
obtain 



m e (T) = (2.484 ± 0.052)T . (5) 

This agrees with the findings of a recent analysis of Polyakov loop correlation func- 
tions, where the electric mass has been determined with the help of the transfer 



matrix approach [16]. We will return to a more detailed discussion of the electric 



mass after having analyzed the magnetic screening mass. 

The magnetic screening mass clearly rises slower than proportional to T. We thus 
may compare the results with the expected behaviour m m ~ g 2 (T)T. In Fig. ?? 
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Figure 3: Electric and magnetic screening masses in units of the temperature versus 
T/T c . The screening masses have been extracted from fits to the gluon correlation 
functions for distances zT > 1. The dashed line shows the fit to the electric mass 
given in Eq.flBI) 



we show the inverse mass, T/m m (T). This should rise logarithmically, if the gauge 
coupling is running according to the leading orders of the renormalization group 
equation, g~ 2 {T) ~ In (T/T c ). We note that this does indeed describe our numer- 
ical results quite well. The spatial string tension is also expected to be sensitive 
to magnetic screening in the high temperature phase and indeed a similar scaling 
behaviour has been observed for it in [|Hf|. 
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Figure 4: Inverse magnetic mass and square root of the spatial string tension (from 
Ref. [17]) versus T/T c . The lines show the fits given in Eqs.flD and (H). 
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We have fitted the magnetic mass with the ansatz 
T 



cg~\T) , (6) 



m m (T) 

assuming the validity of the two-loop formula for g 2 {T) with a free A parameter, 

^ 2 (T) = T ^lnT/A m + i ^ln(21nT/A m ) . (7) 
From a fit of the numerical results for T > 2T C we find 



m m (T) = (0.466 ± 0.015)g 2 (T)T , (8) 

with A m = 0.262(18)T C . We note that our result for the magnetic mass is about 
twice as large as the earlier numerical results |7], ||. This may not be too surprising 
in view of the rather small lattices used in these first studies^. Our result also is 
larger than analytic results based on a semiclassical approximation []nj and gap 
equations || ||. 

It is interesting to relate the magnetic mass to the spatial string tension, which 
also has been found to scale like g 2 T [[HJ. The temperature dependent running 



coupling has been found to be somewhat smaller in the case of v/trj, which results 
from the different constant terms that went into the definition of the A-parameter. 
We may, however, express the result for the spatial string tension in terms of the 
coupling g 2 (T) determined here from the scaling of the magnetic mass. This yields 

yfcT s = (0.368 ± 0.017) {1 — (0.113 ± 0.022)g 2 (T)} g 2 (T)T . (9) 

Combining this with Eq. (§) we find 

m m (T) = (1.27 ±0.10) {1 + (0.113 ± 0.022)# 2 (T)} y/a~ 8 . (10) 

The spatial string tension has also been found to be quantitatively closely related to 

the string tension, y/o^, of the three-dimensional SU(2) gauge theory. We thus can 

b The analysis performed in Ref. could only probe distances zT < 0.5 and the analysis in 
Ref. H has been performed at a value of the gauge coupling which was too large to keep the system 
confined in spatial directions. 
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use the above relation to compare the result for the magnetic mass of the (3+1)- 
dimensional SU(2) gauge theory at finite temperature with the mass gap of the 



three-dimensional SU(2) theory. The latter is found to be about [IS], i.e. 

(4-5)-times larger than the magnetic mass found here. It will be interesting to also 
analyze, in the future, the effective gluon mass in the three-dimensional gauge theory 
within the approach discussed here. This will clarify whether the gluon masses can 
be interpreted as constituent masses of the glueballs and to what extend the thermal 
gluon mass in (3+1) dimensions and the gluon mass in three dimensions are related 
in a similar way as cr s (T) and 03. 

In view of the discussion of the magnetic mass we shall now return to the result for 
the electric mass given in Eq. (|5]). The running coupling extracted from the magnetic 
mass varies by a factor two in the temperature interval studied by us, g 2 (2T c ) ~ 4.1 
and g 2 (16T c ) ~ 2.2 (similar values have been obtained from J~al). From the leading 



order perturbative result, m e {T) = y2/3g(T)T, we thus would expect a 30% drop 
of m e /T in this temperature interval. Moreover, using the perturbative form at 
our largest temperature, we would conclude that the coefficient in front of g{T)T 
comes out to be a factor two larger than suggested by leading order perturbation 
theory. Such an enhancement is also found when higher order effects, including 
the contribution from a non-vanishing mass, are taken into account through re- 
sumed perturbation theory ||. It thus seems that indeed non-perturbative effects 
give large contributions to the electric screening mass and eliminate the leading 
g(T)-dependence. It will be interesting to perform calculations at even higher tem- 
peratures in order to see whether contact with leading order perturbative behaviour 
can be made at all in the electric sector. Work in this direction is in progress. 
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